UNCLASSIFIED 


AD  NUMBER 


ADB075698 


NEW  LIMITATION  CHANGE 
TO 

Approved  for  public  release,  distribution 
unlimited 


FROM 

Distribution  authorized  to  U.S.  Gov't, 
agencies  only;  Test  and  Evaluation;  31  DEC 
1981.  Other  requests  shall  be  referred  to 
Commanding  Officer,  Naval  Research 
Laboratory,  Washington,  DC  20375. 


AUTHORITY 


NRL  [IPO-11A4]  ltr  dtd  8  Aug  1997 


THIS  PAGE  IS  UNCLASSIFIED 


OTIS  FILE  COPY  ADB075698 


XI 


NRL  Report  8568 
(Revised  2/22/83) 


A  Stochastic  Simulation  of  Low  Grazing 
Angle,  Forward  Scatter,  Oyer-Water 
Multipath  Effects 

D.  Y.  Northam 

Advanced  Techniques  Branch 
Tactical  Electronic  Warfare  Division 


December  31,  1981 


Distribution  limited  to  U.S.  Government  Agencies  only,  test  and  evaluation,  December  1981.  Other  requests  for  this  document 
must  be  referred  to  the  Commanding  Officer,  Naval  Research  Laboratory,  Washington,  D.C.  20375. 

83  08. 10  623 


SECURITY  CLASSIFICATION  OF  THIS  PACE  fWi on  D.f*  Entotod) 


REPORT  DOCUMENTATION  PAGE 

HEAD  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

1.  REPORT  NUMBER  2,  COVT  ACCESSION  NO. 

NRL  Report  8568  Ah''$07<S t Cfft 

3.  RECIPIENT'S  CATALOG  NUMBER 

4.  TITLE  (ond  Subtitle) 

A  STOCHASTIC  SIMULATION  OF  LOW  GRAZING 

ANGLE,  FORWARD  SCATTER,  OVER-WATER 

MULTIPATH  EFFECTS 

5.  TYPE  OF  REPORT  i  PERIOD  COVERED 

A  final  report  on  one  phase  of  a 
continuing  NRL  problem. 

6.  PERFORMING  ORG.  REPORT  NUMBER 

7.  author^; 

D.  Y.  Northam 

8.  CONTRACT  OR  GRANT  NUMBERfcJ 

9.  PERFORMING  ORGANIZATION  NAME  ANO  AOORESS 

Naval  Research  Laboratory 

Washington,  DC  20375 

10.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  8  WORK  UNIT  NUMBERS 

57-0481-E-l 

57-0535-0-1 

1 1.  CONTROLLING  OFFICE  NAME  ANO  ADDRESS 

12.  REPORT  DATE 

December  31,  1981 

13.  NUMBER  OF  PAGES 

20 

14.  MONITORING  AGENCY  NAME  »  AOORESSfK  dllloront  ttom  Controlling  Olllco) 

15.  SECURITY  CLASS,  (ot  thh  fpott) 

UNCLASSIFIED 

15a.  DECLASSIFICATION/ DOWNGRADING 
SCHEDULE 

IS.  DISTRIBUTION  STATEMENT  (o  1  this  11. port) 

Distribution  limited  to  U.S.  Government  agencies  only;  test  and  evaluation;  Dec.  1981. 

Other  requests  for  this  document  must  be  referred  to  the  Commanding  Officer,  Naval 

Research  Laboratory,  Washington,  DC  20375. 

17.  DISTRIBUTION  STATEMENT  (o  1  tho  .Stlr.cl  on  to  rod  In  Bloch  JO,  II  dllloront  Irom  Roporl) 

18.  SUPPLEMENTARY  NOTES 

19.  KEY  WORDS  (Continue  on  tovorto  old*  It  nocoooory  and  Identity  by  bloek  numbor} 

Multipath 

Simulation 

Stochastic  models 

20.  ABSTRACT  /Conf/nu#  on  old*  It  noc •»*sry  and  Identity  by  block  nuafbor) 

This  report  describes  a  stochastic  model  and  a  simulation  of  low  grazing  angle,  forward  scatter, 
over-water  multipath  applicable  to  microwave  frequencies.  The  model  is  based  primarily  on  the 
work  of  Charles  Beard,  using  a  coherent  and  an  incoherent  vector  summation  to  describe  the  reflected 
signal.  A  second  order  stochastic  model  and  an  associated  stochastic  simulation  are  described. 
Controlling  parameters  are  sea  state,  radar  frequency,  transmitter/receiver  geometry,  and  radar 
antenna  patterns.  A  model  is  suggested  for  the  incoherent  spectrum  when  one  or  both  radars  are 
moving. 

DO  ^an“3  1473  EDITION  OF  1  NOV  65  IS  OBSOLETE 


j  SECURITY  CLASSIFICATION  OF  THIS  PACE  (Wh *n  D«t«  Bnftod) 

t 

/  i 


CONTENTS 


INTRODUCTION .  1 

MODEL  AND  SIMULATION  CONSTRAINTS  .  1 

MULTIPATH  MODEL .  2 

Model  Summary  .  2 

Coherent  Component  .  3 

Incoherent  Component .  4 

Extrapolation  of  the  Coherent  and  Incoherent  Coefficients  .  5 

Reflection  Coefficient  for  a  Smooth  Plane  Surface .  6 

Receive  Signal  Spectrum  .  8 

Generation  of  the  Incoherent  Processes .  11 

Antenna  Pattern  Effects  .  13 

Moving  Platform  Multipath .  14 

SUMMARY  .  15 

ACKNOWLEDGMENTS .  16 

REFERENCES  .  16 


'Un->  Ko., 

'  l  v 

'-  •»  K 

\ 


n 


n 


m 


A  STOCHASTIC  SIMULATION  OF  LOW  GRAZING 
ANGLE,  FORWARD  SCATTER,  OVER-WATER 
MULTIPATH  EFFECTS 


INTRODUCTION 

One  of  the  major  concerns  of  the  Tactical  Electronic  Warfare  Division  (TEWD)  of  the  Naval 
Research  Laboratory  is  evaluating  the  performance  of  electronic  countermeasures  (ECM)  in  defending 
ships  from  anti-ship  cruise  missile  (ASCM)  attacks.  In  the  Advanced  Techniques  Branch  (Code  5750), 
digital,  pulse-by-pulse  simulations  of  missile-ship  engagements  are  used  to  study  in  detail  the  response 
of  missiles  to  ECM  during  such  engagements.  Because  these  engagements  take  place  at  sea,  over-water 
multipath  is  often  an  important  factor  to  consider  in  such  studies.  For  this  reason,  Code  5750  has 
developed  and  implemented  a  pulse-by-pulse,  digital  simulation  of  low  grazing  angle,  over-water,  for¬ 
ward  scatter  multipath  applicable  to  microwave  radar  frequencies. 

MODEL  AND  SIMULATION  CONSTRAINTS 

The  over-water  multipath  phenomenon  is  of  course  an  extremeiy  complex  one  to  describe  in 
detail.  Therefore,  it  was  important  to  isolate  the  dominant  effects  of  the  phenomenon  on  the  particular 
ECM  studies  that  Code  5750  anticipated  conducting.  Further,  it  was  important  to  develop  mathematical 
models  that  are  amenable  to  validation  (by  validation  we  mean  that  the  model  accurately  describes  the 
effects  to  second-order  statistics).  Finally,  the  model  had  to  lead  to  a  digital  simulation  that  would  be 
used  as  a  submodule  in  larger,  more  comprehensive  simulations  whose  purpose  would  be  ECM  evalua¬ 
tion.  These  general  requirements  lead  to  more  specific  requirements  or  constraints.  Because  the  simu¬ 
lation  was  to  be  a  part  of  several  larger  pulse-by-pulse  simulations,  it  had  to  be 


(1)  able  to  generate  multipath  effects  as  a  time  series  (pulse-by-pulse), 

(2)  fast  running,  and 

(3)  modular,  for  general  applicability  and  clarity. 

Because  the  applications  had  to  deal  with  understanding  detailed  aspects  of  radar  systems  performance, 
the  simulation  had  to 

(4)  generate  the  complex  (amplitude  and  phase)  multipath  effects,  and 

(5)  generate  receive  signals  (affected  by  multipath)  accurate  to  at  least  second  order  statistics 
(i.e.,  density  function,  location  and  spread  parameters,  and  power  spectral  density). 

Because  the  mathematical  model  had  to  be  validated,  its  parameters  needed  to  be 

(6)  as  few  as  possible,  and 

(7)  experimentally  measurable. 
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These  constraints  have  led  to  the  development  of  a  second-order  stochastic,  phenomenological  model, 
that  is,  a  model  that  is  based  on  physically  measurable  parameters  and  that  describes  the  extremely 
complex  physical  processes  (the  generation  of  incoherent  scattering  by  a  rough  sea  surface)  stochasti¬ 
cally. 

MULTIPATH  MODEL 

The  stochastic  model  which  we  have  developed  is  based  primarily  on  the  work  of  Beard.  Katz,  and 
Spetner  [l-'j.  John  Baras  developed  the  original  version  of  the  stochastic  model  described  in  this 
report  and  has  also  developed  several,  more  sophisticated  (detailed)  multipath  effects  models  (4)  than 
the  one  presented  in  this  report.  Simulations  of  the  more  sophisticated  models  have  not  been 
developed  because  for  each  of  our  application  studies  to  date,  we  have  judged  the  Beard  model  to  be 
adequate  in  its  description  of  multipath  effects. 


Model  Summary 

The  simulation  model  is  based  on  the  concept  of  vectors  representing  the  electric  field  strength  at 
the  receive  antenna.  This  field  strength  is  modeled  as  the  sum  of  two  components,  a  direct  path  vector 
and  a  reflected  vector. 

Figure  1  illustrates  and  defines  the  individual  vector  components  of  the  model.  The  direct  signal 
D  (that  received  in  the  absence  of  multipath  effects)  is  taken  as  the  referen£e  vector.  The  reflected  sig¬ 
nal,  R,  is  the  sum  of  two  other  vectors,  a  coherent  (relative  to  D)  vector,  C,  and  an  incoherent  vector, 
/.*  The  coherent  component  is  the  specularly  reflected  vector  which  would  be  received  if  the  sea  was 
perfectly  smooth,  modified  by  a  scattering  coefficient.  The  incoherent  component  accounts  for  the  time 
varying  effects  of  the  sea  on  the  receive  signal  (in  particular,  the  spectrum).  The  non-time-varying  pro¬ 
perties  of  the  reflected  components  are  primarily  functions  of 

(1)  the  electrical  properties  of  the  sea, 

(2)  the  geometry,  and 

(3)  the  radar  wavelength,  polarization,  and  antenna  patterns. 

The  total  (vector  sum)  receive  signal  is  represented  by  the  vector  T. 


Fig.  1  —  Receive  signal  vector  summation 


Consider  transmit  and  receive  radar  antennas  fixed  relative  to  the  earth,  a  distance  R  apart,  and  at 
heights  above  the  mean  sea  level  of  hT  and  /i^,  respectively  (Fig.  2).  If  the  geometry  does  not  change 

The  coherent  and  incoherent  vectors  are  also  referred  to  as  the  specular  and  diffuse  sectors,  respectively. 
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with  time  (an  assumption  that  we  make  to  develop  the  fundamental  model)  then  the  only  significant 
variation  in  the  receive  signal  is  due  to  sea  surface  motion.  This  motion  directly  controls  the  receive 
signal  variation,  and  we  use  a  characterization  of  its  second-order  properties  to  determine  the  receive 
signal’s  second-order  properties.  For  our  purposes,  the  sea  surface  motion  can  be  adequately  character¬ 
ized  by  a  simple  representation  of  the  wave  height  spectrum  and  a  surface  roughness  parameter  g. 
These  two  observables  are  the  most  fundamental  parameters  of  the  model.  The  roughness  parameter 
used  is  defined  by 


/?  sin 
X 


(1) 


where 


h  --  RMS  wave  height  (standard  deviation) 


i ji  =  grazing  angle 

X  =  radar  wavelength. 


The  following  sections  describe  in  detail  the  various  components  of  the  simulation  model. 


Coherent  Component 

The  coherent  vector  is  described  by 

,2nbR 

C~\D\rPce  x  , 

where 


D  —  direct  signal 

T  =  smooth  sea  reflection  coefficient  (complex) 

pc  ~  coherent  scattering  coefficient 

A R  =  direct  and  reflected,  path  length  difference 

X  =  radar  wavelength. 


(2) 
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Antenna  pattern  effects,  not  shown  in  Eq.  (2),  are  discussed  in  the  section  devoted  to  this,  and  a  sum¬ 
mary  of  the  calculation  of  T  is  also  given  in  a  later  section.  Figure  3  is  a  plot  of  three  models  that  have 
been  developed  to  describe  pc.  Ament’s  model  [5]  is 

Pe  =  e"2<2**>\  0  <  g  <  0.3.  (3) 


Fig.  3  —  Three  models  for  the  coherent  scattering  coefficient  (from  (61). 
Solid  line  is  theory  (Brown  and  Miller),  dashed  line  is  approximation  to  ex¬ 
perimental  curve  (Beard),  and  dotted  line  is  theory  (Ament). 


Brown  and  Miller  [6]  have  derived  an  expression,  using  geometrical  optics  which  is  a  better  fit  to 
Beard’s  data: 

Pc  -  e-™”s)2I0(2(2Trg)2),  0  <  s  <  0.3,  (4) 

where  /00c)  is  the  modified  Bessel  function  Jo(ix).  The  third  model  is 

expl— 2(27t£)2],  0  <  g  <  0.1 

Pc=  0-812537  ’  (5) 

1  +  2(2tts)2'  0.1<g<0.3 

which  is  our  approximation  to  the  empirical  fit  made  by  Beard  to  experimental  data  [31.  Models  (4) 
and  (5)  are  implemented  in  the  simulation,  and  the  user  has  the  option  of  specifying  which  model  to 
use  (when  he  initializes  the  simulation). 

Incoherent  Component 

The  incoherent  vector  is  assumed  to  be  the  resultant  sum  of  vectors  from  many,  independent, 
random  scatterers.  Therefore  it  is  represented  as  the  sum  of  two  zero  mean  Gaussian  processes,  Ip  and 
Iq  which  are  orthogonal  (Fig.  1)^  As  in_Ref.  1,  we  assume  that  IP  and  Iq  are  independent  although  it 
is  known  [7]  that  coi  relation  of  Ip  and  Iq  can  occur  due  to  wave  front  sphericity.  The  validity  of  this 
assumption  is  further  strengthened  by  Ref.  8  where  it  was  found  thatjor  data  collected  at  the  Gulf  of 
Mexico,  lP  and  Iq  were  independent.  For  convenien£e,  the_phase  of  IP  is  assumed  to  be  that  of  C.  It 
is  further  assumed  that  the  power  spectral  density  of  IP  and  Iq  are  identical  since  we  know  of  no  physi¬ 
cal  reason  to  assume  otherwise.  Therefore  we  have 

7-  |5|r  (7p  +  7q).  (6) 
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The  standard  deviation  of  the  processes  are  identical  and  are  determined  by  the  incoherent  scattering 
coefficient,  p;.  Figure  4  shows  experimental  data  collected  and  a  curve  fit  to  that  data  by  Beard  [3] 
describing  p/  as  a  function  of  the  roughness  parameter. 


Fig.  4  —  Incoherem  scattering  coefficient  or  normalized  standard  deviation  (from  (31) 


In  the  simulation,  the  expression 

7  =  \d\  \r\S2p,ux  +  ji2)  (7) 

is  used  to  generate  1.  I\  and  /2  are  independent  0,  1)  processes  (normalized  Gaussian).  The  V2 
term  arises  because  the  sum  of  two  N(0,  p ,)  processes  is  a  V2p//V(0,  1)  process.  The  receive  signal 
spectrum  will  be  discussed  later  in  detail. 

Extrapolation  of  the  Coherent  and  Incoherent  Coefficients 

The  above  models  of  pc  and  p,  are  validated  for  g  <  0.3  [3]  and  at  present,  the  author  is  unaware 
of  any  validated  models  for  g  >  0.3.  This  raises  two  questions:  first,  what  are  good  models  of  pt  and 
Pz  for  g  >  0.3  and  second,  given  the  present  models,  how  should  one  simulate  a  transition  from  g  < 
0.3  to  g  >  0.3.  When  a  validated  model  becomes  available  it  can  be  directly  implemented  in  the  simu¬ 
lation.  For  our  applications,  the  second  problem  occurs  only  as  a  result  of  the  grazing  angle  increasing 
during  a  simulation  (\  and  hRMS  are  fixed  during  the  simulation).  This  increase  in  g  is  usually  also 
accompanied  by  antenna  gain  reductions  in  the  reflected  signal  direction  which  further  reduce  the  mul¬ 
tipath  effects.  Therefore  for  our  purposes,  we  assume  that  whatever  pc  and  p/  are  when  g  increases 
beyond  0.3.,  they  are,  for  practical  purposes,  nonincreasing  and  small.  In  the  simulation  we  have  extra¬ 
polated  the  pf  and  p/  models  as  described  in  the  following  paragraph. 

Assuming  minimal  multipath  effects  for  g  >  0.3,  we  suggest  three  practical  approaches  to  simulat¬ 
ing  the  transition  of  pc  and  p/  through  g  =  0.3.  First,  one  can  simply  not  calculate  multipath  terms  for 
g  >  0.3.  Because  this  introduces  discontinuities  in  the  pulse-by-pulse  signal  synthesis  (at  g  =  0.3)  and 
because  there  are  probably  some  "residual''  multipath  effects  for  large  values  of  g ,  we  have  taken  a 
second  approach.  Currently  in  the  simulation  we  make  smooth  extrapolations  of  both  pt  and  py  pt  is 
extrapolated  by  extending  (4)  and  (5)  to  all  values  of  g ,  and  p;  is  a  linear  extrapolation  (Fig  5)  of  the 
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linear  fit  in  Fig,  4,  The  third  method  is  to  smoothly  extrapolate  the  two  coefficients  to  zero  at  some  /; 
“  A' max  >  0,3  and  not  calculate  any  multipath  terms  (which  would  be  zero)  for  #  > 

Reflection  Coefficient  for  a  Smooth  Plane  Surface 

Here  we  summarize  the  model  of  the  reflection  coefficient  following  Chapters  2  and  5  of  Kef.  9, 
There  arc  two  expressions  for  the  coefficient  depending  on  the  radar  polarization: 


sin  ih  -  y/T^-co^ili 


sin »//  +  V«(-  -  cos2  ih 


for  vertical  polarization  and 


•II  /— - a — 

sin  ih  +  v«,  -  cos'  ih 

for  horizontal  polarization  where  »//  is  the  grazing  angle  and  t,  is  the  complex  dielectric  constant  which 
we  approximate  by 

«,  -  t  -  60\irl  (10) 


where 


<  is  sea  dielectric  constant, 

A  is  radar  wavelength  (m),  and 

<t  is  sea  conductivity  (mho/m). 

Based  on  Table  5,1  of  Ref,  9,  we  have  chosen  to  use  seawater  values  of  a  —  4,3  mho/ m  and  «  ■*  80, 
Simulations  of  (8)  and  (9)  have  shown  that  for  our  applications,  variations  in  1" //  and  IV  Cue  to  varia¬ 
tions  in  t  and  <r  about  these  values  arc  not  significant.  Figures  6  through  9  illustrate  the  character  of 
the  amplitude  and  phase  of  IV  and  V n  as  a  function  of  grazing  angle  for  fixed  A,  <r,  and «. 
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Fig.  9  —  Angle  of  l'i ;  X  =■  2.5  cm,  it  *■  4.3  mho/m,  e  ■»  80 


The  calculation  of  T  is,  relatively,  very  time  consuming  and  variations  in  T  due  to  small  changes 
in  grazing  angle  are  also  small.  Therefore  we  include  in  the  simulation  a  test  to  decide  when  to  recalcu¬ 
late  F  rather  than  performing  a  pulse-by-pulse  calculation.  This  is  particularly  efficient  for  stationary 
platform  studies.  The  thresholds  (indicating  significant  grazing  angle  change)  used  were  determined  by 
examining  F  where  the  rate  of  change,  of  F  was  most  significant.  A  threshold  of  1°  change  was  chosen 
for  l" n  (somewhat  arbitrarily  due  to  the  minimal  variations  in  the  amplitude  and  phase  of  F,,).  An 
examination  of  the  T  v  angle  curve  indicates  that  as  much  as  a  180°  phase  change  can  occur  for  about  a 
1°  grazing  angle  change.  Therefore,  a  value  of  0.01°  is  used  in  the  simulation  as  an  update  threshold 
for  calculating  T,/ 

Receive  Signal  Spectrum 

Up  to  this  point  we  have  described  only  the  first-order  properties  of  the  receive  signal  (RS).  To 
complete  the  second-order  model,  we  need  a  description  of  the  RS  spectrum  (or  autocorrelation  func¬ 
tion). 


For  our  applications,  we  place  two  requirements  upon  the  character  of  the  RS  variations: 

(1)  for  fixed  geometry  (Fig.  2),  the  RS  must  be  statistically  stationary,  and 

(2)  wave  height  directionality  must  not  be  a  parameter  of  the  RS  variations  (our  applications 
require  an  "average  effect"  with  respect  to  wave  directionality). 

These  requirements  will  be  assumed  throughout  the  following  paragraphs. 

When  the  geometry  is  fixed,  the  fluctuations  in  the  RS  are  due  to  the  motion  of  the  sea  surface 
Therefore,  the  various  characteristics  of  the  RS  fluctuations  should  be  directly  related  to  the  associated 
characteristics  of  the  sea  surface  variations.  This  intuitive  idea  has  been  quantitatively  verified 
[2,8,10,11].  In  particular,  the  experimental  work  described  in  Ref.  10  indicates  that  the  RS  spectrum 
is,  in  fact,  .rongly  related  to  the  wave  height  (WH)  spectrum.*  In  particular,  both  spectra  have  been 
shown  to  possess  at  least  for  relatively  low  sea  stales,  the  same  qualitative  shape.  Each  spectrum  has 
what  may  roughly  be  described  as  a  "bandpass  shape",  that  is  they  can  be  characterized,  in  large  meas¬ 
ure,  by  specifying  a  center  frequency,  a  bandwidth,  and  a  process  variance. 

*A  detailed  quantitative  analysis  of  these  data  has  not  been  conducted  to  date.  However  it  is  anticipated  that  such  work  will  be 
undertaken  in  the  near  future 
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For  our  simulation  we  have  chosen  to  generate  the  spectrum  of  the  RS  process  by  modeling  the 
spectra  of  the  incoherent  components  as  bandpass  spectra.  This  was  done  for  two  reasons.  First, 
bandpass  processes  can  be  easily  and  efficiently  generated  digitally  and  second,  although  it  is  known  that 
the  WH  and  RS  spectra  are  related,  a  precise  quantitative  description  of  the  RS  spectrum  is  not  yet 
known.  However,  the  following  approximate  results  have  been  found  experimentally: 

Ccrkt.r 

(1)  the  upper  3  dB  frequency  of  the  RS  spectrum  is  linearly  related  (Fig.  10)  to  the 
M  frequency  of  the  WH  spectrum  [2],  and 

(2)  the  RS  spectrum  center  frequency  is  approximately  equal  to  the  WH  spectrum  center  fre¬ 
quency  [10]. 


0  0.2  0.4  0.6 

g 


Further,  by  using  a  series  expansion  of  the  magnitude  of  the  total  receive  signal  T  (Fig.  1),  Beard  and 
Katz  [2]  have  shown  that 

(3)  when  D  and  C  are  in  phase,  to  first  order,  the  spectrum  of  \lP\  is  equal  to  the  spectrum  of 

171. 

Our  approach  to  generating  the  spectrum  of  the  RS  process,  based  on  the  above  three  results,  is. 

(1)  Choose  a  sea  state. 

(2)  Based  on  the  chosen  sea  state,  pick  a  representative  WH  spectrum. 

(3)  Calculate  a  center  frequency,  fw,  for  the  WH 

spectrum. 

(4)  To  generate  the  RS  process  with  appropriate  spectrum,  generate  IP  and  In  processes 
which  have  bandpass  spectra  each  with  center  frequency  f\V  (from  result  (1)  above), 
upper  3  dB  frequency  fRS( 0.5)  (from  result  (2)  above),  and  variance  %/2  •  p,  (from  the 
section  describing  the  incoherent  component). 

Summing  D  and  C  vectors  with  lP  and  Iq  as  described  above  produces  a  power  process  !  T P  whose 
second  order  statistics  match  closely  the  second-order  statistics  of'experimental  data  collected  by  Code 
5750. 
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Although  various  models  of  WH  spectra  exist,  for  our  purposes  we  need  only  a  model  that  will 
yield  an  estimate  of  fw  for  a  given  sea  state  or  a  given  wind  speed.  For  this  we  use  the  Neumann  spec¬ 
trum  (12]  which  is  defined  by 


iKf) 


£1 


(ID 


where 


V 

f 

C> 

C2 


Se 


=  wind  speed  (knots) 


=  frequency  (Hz) 
=  0.779 


=  32.2  ft/s2. 


The  general  form  of  this  spectrum  is  shown  in  Fig.  11  Equation  (11)  can  be  used  to  obtain  the  center 
frequency  fw  by  setting  the  derivative  of  A^(J)  to  zero  and  solving  for/J(/.  This  yields 


fw  ~ 


2.476 
V  ' 


(12) 


The  variance  of  the  WH 'process  (required  for  the  calculation  of  g)  can  be  estimated  by  integrating  the 
assumed  (Neumann)  WH  spectrum.  Using  the  standard  symbol  E ,  the  variance  is  [13] 


where 


and 


\ 


E 


0.242 


> 


V  =  wind  speed  (knots).  / 


(13) 
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From  Eq.  (13)  we  obtain  the  RMS  wave  height  (WH  standard  deviation)  as 

hRXIS  =  0.00110  V*2.  (14) 

Finally,  we  use  the  power-point/center-frequency/sea-roughness  curve  of  Fig.  12  to  obtain  fRS( 0.5). 
The  straight  line  portion  of  this  curve  is  based  on  Fig.  10  ( g  ^  0.1).  We  have  implemented  a  smooth 
interpolation  (using  a  C,  +  C2g  3  fit)  from  5  =  0.1  to  g  —  0.0.  This  allows  continuous  variations  of 
fRS( 0.5)  in  the  simulation. 


VO 

2  '  10 


0.4  0.5 


Fig.  12  —  Calculation  of/RS  (0.5)  from  sea  roughness  and  wave 
spectrum  center  frequency 


In  general,  our  studies  use  one  of  three  nominal  sea  conditions:  slight,  moderate,  or  rough.  To 
obtain  /.’ra/s.  /r$( 0.5),  and  fn  for  a  given  sea  state  we  first  chose  a  wind  speed  appropriate  to  the  sea 
state,  then  calculate  the  above  parameters.  The  values  of  wind  speed  that  we  typically  use  along  with 
values  of  associated  parameters  are  shown  in  Table  1. 

Table  1  —  Typical  Values  of  Sea  Surface  Parameters  for 
Three  Sea  Conditions 


Sea  state 


Wind  speed 
(knots) 


BWw 

fw 

(Hz) 

(Hz) 

Slight  (SS2) 
Moderate  (SS3) 
Rough  (SS4) 


Generation  of  the  Incoherent  Processes 


The_  fundamental  stochastic  processes  that  we  generate  in  the  simulation  are  the  incoherent  com¬ 
ponents  IP  and  Iq.  As  stated  before,  for  the  fixed  platform  model  the  processes  will  be  wide-sense  sta¬ 
tionary,  zero-mean  Gaussian,  and  independent.  Because  the  processes  have  bandpass  spectra,  they  will 
have  autocorrelation  functions  of  the  form 


R  (t)  =  exp  (-ar)  cos  tor. 
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This  form  and  the  requirement  that  a  fast  digital  simulation  be  developed  led  us  to  use  first-order, 
Gauss-Ma.kov  random  sequences  to  generate  the  processes.  A  detailed  analysis  of  such  sequences  is 
given  in  Ref.  14. 

We  begin  the  description  of  the  process  generation  by  describing  the  generation  of  an  exponen¬ 
tially  correlated,  Gaussian  sequence  XK  (see  also  Ref.  14).  Let  WK  be  a  Gaussian,  white  noise 
sequence  with 


£{^1  =  0 

,  ,  k  =  1 
e{wkw,)=  d>k9il 


Then  the  values  of  the  process  at  time  K  can  be  approximated  by 


XK  -  AXk.x  +  BWk 

where  ,  „ 

£  W  =  0  '  (17> 

£(T02)  =  Rx  (0)  =  R0 

We  need  to  find  values  for  A  and  B.  By  calculating  the  covariance  of  XK,  one  can  show  [14]  that  the 
covariance  will  be  stationary  if 

0  <  A  <  1,  and 

d2  •  (18) 

R‘-jh 

This  defines  B  in  terms  of  /?0  and  A.  Next,  we  require  that 

Rq  exp  (-2Ar)  (19) 

where  Ar  is  the  simulation  step  size.  We  also  have 

£{,%,}  =  E{X0(AX  o  +  BW A()} 

=  AE{X%)  +  BE[XoW^ 


so  that  A  is  determined: 


R0  exp  (-aAf)  =  AR0, 


A  =  exp  (-aAt). 

Finally,  using  £0  from  Eq.  (17)  the  XK  process  is  stationary  since 


and  [14], 


£{**}  =  0, 


E{XkXkJ  =  yfyy  |a|'. 


We  now  have  a  method  of  generating  a  first-order,  Gauss-Markov  sequence  with  autocorrelation 
Rq  exp  (-a  /),  using  a  step  size  A  / : 
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where 


Xk  =  AXk_,  +  BWk 
A  =  exp  (— 2A  t) 

B  =  V*o0  "  ^2) 

£{*0}  =  0 

£U02}  =  £o- 


(22) 


To  complete  the  incoherent  process  generation,  we  need  to  transform  XK  to  a  process  with  a 
bandpass  spectrum.  Let 

Yx  —  X i  x  cos  oit  4-  X2  k  sin  0 it,  (23) 

where  XXK  and  X2K  are  independent  and  as  above  (i.e.,  N(0,  y/R~0),  stationary,  with  RX[(r)  —  Rx2 (T) 
=  Rx(t))  and  t  —  KRt.  The  YK  process  is  also  Gaussian  because  it  is  a  memoryless  linear  transforma¬ 
tion  of  Gaussian  processe:.  Further,  it  is  stationary  with  the  desired  autocorrelation  because: 

Ely*)  —  cos  o)tE[X\  *}  4-  sin  cotE{X2  %) 

=  0, 


and 

E{Yf;YK+,}  =  E{{Xxk  cos  ait  4-  X2K  sin  cot) 

(Xi,K+i  cos  a>(r  +  r)  +  X2K+i  sin  w(r  +  r))) 

=  Rx  (r )  cos  Oi  t  COS  0)  (t  4-  t) 

+£{Z1  ,^^2,^+/}  cos ait  sin«0  +  t) 

+£{^2X^1.^+;}  sin  cor  cose o(r  4-  r) 

+Rx(t)  sin uit  sin  co(r  4-  r) 

=  Rx{t)(cos  ait  cos  c o(r  +  t)  +  sin  cor  sin  co(r  +  r)) 

=  RX(t)  COS  OIT.  (24) 

We  now  have  a  generation  procedure  for  the  two  incoherent  processes.  Using  four  XK  processes  we 
generate 

Ipx  =  Pi[X\  X  cos  (cotfAr)  4-  X2x  sin  (&> /if A r ) 3 
Iq.k  “  pAX^k  cos  (c oKkt )  +  X4K  sin  (co/TAr)] 

where  we  have  set  £0  ~  L 


Antenna  Pattern  Effects 


Our  ECM  applications  simulations  allow  the  use  of  various  antenna  patterns  depending  on  the 
particular  radar  being  simulated.  Because  of  this,  the  multipath-effects  simulation  does  not  include 
antenna  pattern  effects  explicitly.  (The  antenna  pattern  and  multipath  simulations  are  separate 
modules.)  These  effects  must  be  introduced  by  the  user  once  he  has  chosen  an  antenna  pattern  and 
determined  a  need  to  include  multipath  effects. 

The  effects  of  antenna  patterns  on  the  received  signal  are  not  well  understood  quantitatively  and 
the  author  is  not  aware  of  any  validated  model  for  these  effects.  However,  a  simple  model  to  account 
for  these  effects  (which  we  have  used  in  the  simulation  to  date)  is  described  in  the  following.  We 
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assume  that  the  two  receive  signals  ( D  and  R)  are  signals  transmitted  from  two  point  sources  located  at 
the  transmit  radar  and  the  sea  surface  specular  point  and  that  they  are  summed  by  the  receive  antenna 
with  the  associated  antenna  gains  (Fig.  13).  Because  both  D  and  R  are  representations  of  electric  fields, 
to  correctly  account  for  antenna  effects  in  such  a  model,  amplitude  gains  must  be  used.  This  is  done  by 
producing  a  complex  multipath  factor  Fmp  normalized  to  the  direct  signal  D  =  D  exp  (J0D ).  It  is 

Fmp  =  F  exp  ( jO )  =  Fr  +  jF,.  (26) 


Assuming  unity  antenna  gains,  the  total  received  electric  field  becomes 

T  =  D  [(1  +  Fr)  +  JF,]  exp  (J0D).  (27) 

Now  suppose  the  transmit  and  receive  antennas  have  direct  path  amplitude  gains  GTD,  and  GRl), 
respectively,  and  reflected  path  amplitude  gains  of  GTR  and  GR  R ,  respectively  (all  gains  may  be  com¬ 
plex).  Then  the  received  electric  field  including  antenna  effects  is  represented  by 

7 'ant  -  DIGr.dGt.d  +  GRiRGT,RFmp]  exp  (J0D ). 

This  simple  model  accounts  for  the  qualitative  fact  that  as  the  antennas  illuminate  less  and  less  of  the 
sea  surface,  the  effect  of  multipath  on  the  received  signal  should  be  reduced.  However,  it  has  been 
shown  by  Beard  that  this  model  is  not  quantitatively  valid  when  the  antennas  do  not  fully  illuminate  the 
first  Fresnel  zone  [3,15],  Even  so,  we  have  used  this  model  of  antenna  patterns  effects  in  the  simula¬ 
tion  for  three  reasons:  at  present  there  exists  no  validated  model  for  these  effects;  this  model  does  cap¬ 
ture  the  qualitative  effect  of  reducing  the  surface  illumination;  and  it  is  a  simple  model  to  simulate 
When  an  improved  model  is  available,  it  will  be  incorporated  in  the  simulation. 

Moving  Platform  Multipath 

The  model  as  described  to  this  point  has  assumed  that  the  transmit  and  receive  antennas  were 
fixed  relative  to  the  earth.*  However,  since  the  situation  where  either  one  or  both  of  the  radars  are 
moving  is  of  prime  interest  to  us  in  our  studies,  we  now  consider  relaxing  the  fixed  antenna  constraint 
The  deterministic  parts  of  the  fixed  model  describe  the  deterministic  moving  platform  multipath 
(MPM)  effects  by  varying  the  geometry  (antenna  locations)  with  time.  This  in  turn  produces  time  vari¬ 
ations  in  the  associated  deterministic  parameters  of  the  MPM  model.  However,  the  random  com¬ 
ponents  of  the  fixed  model  (the  incoherent  processes)  are  controlled  by  parameters  that  are  determined 
by  reflections  from  a  fixed  area  of  the  sea  surface.  If  either  or  both  of  the  platforms  move,  then  the 
area  of  reflection  will  no  longer  be  fixed. 

Because  the  author  is  unaware  of  any  validated  model  for  the  incoherent  component  for  MPM, 
we  have  developed  a  heuristic  approach  to  modeling  this  component,  which  we  will  now  discuss  First, 
it  appears  reasonable  to  assume  that  lP  and  Iq  will  remain  independent,  N{0,pt)  processes  and  that  p  / 


•Relative  motion  has  been  anticipated  in  the  previous  sections;  e.g.  the  p;  and  pc  extrapolations. 
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will  remain  a  function  of  only  the  sea  roughness  factor.  The  incoherent  spectrum  is  directly  deter¬ 
mined  by  the  sea  surface  structure.  Because  this  structure  is  bandpass  in  nature  (wavelengths  occur 
within  certain  ranges  for  a  given  sea  state)  and  because  we  assume  that  the  sjrface  is  statistically  sta¬ 
tionary,  we  expect  the  MPM  incoherent  spectrum  to  remain  bandpass  in  shape.  The  major  area  of  sea 
surface  contributing  to  the  multipath  signal  is  the  area  within  the  first  several  Fresnel  zones  [3,8]. 
Therefore,  if  the  platforms  are  moving  relative  to  one  another,  these  zones  are  also  changing  location 
and/or  size  proportional  to  the  rate  of  change  of  geometry  of  the  platforms.  If  the  zones  are  moving 
rapidly  relative  to  the  sea  surface  motion,  then  we  can  view  the  sea  as  fixed  for  a  given  period  of  time. 
If  during  this  period  the  zones  (sea  surface  area)  change  enough  to  cause  a  decorrelation  of  the 
received  incoherent  component,  then  one  can  associate  the  spectral  bandwidths  of  the  incoherent 
processes  with  a  decorrelation  distance  (which  translates  to  a  decorrelation  time)  due  to  platform 
motions.  Viewing  the  sea  surface  as  fixed,  with  periodic  structure  centered  about  some  dominant  sea 
wavelength  (a  "typical"  crest-to-crest  length),  we  further  suggest  that  the  incoherent  spectrum  will  have 
a  dominant  (center)  frequency  inversely  proportional  to  the  time  required  for  the  specular  point  to 
move  from  "dominant  wave  peak"  to  "dominant  wave  peak."  Finally,  we  assume  that  the  energy  of  the 
incoherent  components  will  remain  approximately  constant  for  constant  sea  roughness  factor. 

Because  our  studies  to  date  involving  MPM  (using,  typically  about  Mach  1  closing  velocities)  have 
not  required  detailed  knowledge  of  the  incoherent  spectrum,  only  the  coherent  induced  nulling,  we 
have  implemented  a  simple  approximation  to  this  spectrum  based  on  the  above  discussions.  Our 
approximation  is: 

(1)  assume  a  bandpass  spectral  shape, 

(2)  calculate  a  center  frequency  proportional  (o  the  ratio  of  specular  point  velocity  and  some 
"typical"  sea  wave  length,  and 

(3)  calculate  a  bandwidth  proportional  to  an  estimate  of  the  Fresnel  zone  decorrelation  time. 

Code  5750  is  investigating  the  possibility  of  conducting  an  experiment  involving  a  fixed  radar  and 
an  airborne  radar  to  collect  data  appropriate  for  further  investigation  of  the  MPM  incoherent  spectrum. 

SUMMARY 

In  this  report  we  have  presented  a  simulation  model  of  the  effects  of  multipath  in  low  grazing 
angle,  over-water,  radar  scenarios.  The  purpose  of  the  simulation  is  to  generate  a  time  series  of  data 
accurate  to  second-order  statistics  in  its  representation  of  experimentally  measured  effects,  in  particular, 
those  effects  measured  by  Code  5750  [10]  which  further  confirm  the  fundamental  vector  representation 
made  by  Beard  et  al.  [1].  The  simulation  model  is  based  on  the  vector  model  of  Beard  et  al.,  and  util¬ 
izes  various  experimental  results. 

There  are  three  major  weaknesses  in  the  model: 

(1)  The  inability  to  accurately  account  for  antenna  pattern  effects  when  the  first  Fresnel  zone 
is  not  uniformly  illuminated. 

(2)  The  extrapolation  of  the  scattering  coefficient  (coherent  and  incoherent)  to  regimes  where 
they  have  not  been  experimentally  verified. 

(3)  The  use  of  a  model  for  the  diffuse  spectrum  which  has  not  been  validated  for  moving 
platforms. 

It  should  be  noted  that  these  weaknesses  can  all  occur  in  moving  platform  scenarios. 
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This  report  describes  an  application  of  ongoing  experimental  and  analytic  research  being  con¬ 
ducted  within  Code  5750.  Such  research  and  changing  application  requirements  will  undoubtedly  lead 
to  modifications  to  the  simulation,  and  so  the  simulation  source  code  has  been  modularly  designed  to 
allow  easy  alterations  of  each  component  of  the  model.  In  particular,  further  analysis  of  Code  5750 
multipath  data  could  yield  an  improved  model  of  the  receive  signal  spectrum.  As  appropriate,  any  such 
modifications  will  be  reported.  Another  report  is  being  prepared  that  contains  the  FORTRAN  source 
code  listings  used  by  Code  5750  to  implement  the  simulation  described  in  this  report. 
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